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I show for the specific case of the scalar field spectrum on regular tessel- 
lations of the sphere that the first two terms of the heat-kernel expan- 
sion are related to the first two terms of the Ehrhart (quasi)polynomial. 
In trying to make this relation precise, I consider degeneracies as parti- 
tion denumerants and show the connection of the group theory expres- 
sions with Popoviciu's theorem and with the notion of Sylvester waves. 
General denumerants are considered and the first wave, W 1 , i.e. the 
polynomial part, is written using the A-genus multiplicative sequence. 
It is pointed out that Sylvester in effect did the same thing and that he 
had also obtained Ehrhart reciprocity. I derive an algebraically neat 
form for the second wave, W 2l which involves the combination of two 
multiplicative sequences. 
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1. Introduction 



Connections between lattice and spectral problems go back over a hundred 
years. I mention only the Weyl conjecture and its physical antecedents (see, for 
example, Baltes and Hilf, [1]). 

Two results, which are curiously similar, are the Ehrhart polynomial and the 
short-time expansion of the heat-kernel. Actually, rather than the heat-kernel, it 
is more appropriate to consider the asymptotic behaviour of a smoothed version, 
iV(A), of the exact eigenvalue counting function, N(X) = J2 k x <x l. 2 As a power 
series in A, this has coefficients simply related to those of the heat-kernel expansion. 
For simplicity of exposition, I restrict initially to two dimensions and give the first 
two terms (which are all I am interested in) explicitly, 

N(u) ~ ^(\M\u 2 ± \dM\u + c + ...^J . (1) 

I have set A equal to u 2 and redefined N(ui) = N(X). 

These first two terms are those relevant for the Weyl conjecture, which states 
that (1) holds for the exact counting function. While this may not be true (actually 
it isn't for the manifolds I look at), (1) does hold for the smoothed version. 

Turning to the lattice side, Ehrhart proved that the number of integer lattice 
points (i.e. elements of Z 2 ) in, and on, a rational polygon, i.e. a polygon whose 
vertices are rational numbers, that has been uniformly dilated by an integer, /, is a 
quasipolynomial in /. 

More precisely let V be a rational polygon 3 and V its closure (i.e. including 
its boundary) and let L(P, /) = %(IV H Z 2 ) be the number of integer points in the 
closure of the dilated polygon IV = {(lx, ly) : (x, y) G V ; / G N}. Then 

L(V,l) = \V\l 2 + Cl (l)l + c (t) (2) 

where c l7 c are generally periodic functions in / (see, for example, Beck and Robins, 
[2] and Wright [3]. 

In the special case that the polygon is an integer polygon the coefficients c x , c 
are constant and, moreover, 

L(VJ) = \V\l 2 + ^\dV\l + c (3) 

2 The eigenvalues, A fe are ordered linearly, with k a counting label. 

3 1 sometimes refer to a disc with polygonal boundary (a 2-polytope) as a polygon, for short. 
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which bears a remarkable similarity to (1). 

The two expressions cannot be related in all circumstances because not all lat- 
tice counting problems have spectral associations. This paper, however, is concerned 
with a very particular case in which they can be more precisely connected. 

In sections 2 and 3, I set up the situation and present an experimental result. 
The later sections are intended to be more exact and give some explanation of the 
result. This paper thus proceeds from the particular to the general. 

2. The factored sphere 

The situation I refer to is very well known and so I need not spend time on 
the underlying mathematical details. They have been dealt with by, for example, 
Gromes, [4], Berard and Besson, [5], Briining and Heintze, [6] and Chang and 
Dowker, [7]. The manifold, M is an orbifold factor of the two-sphere, S 2 /r, where 
T is the reflective symmetry group of one of the regular solids, a finite subgroup of 
0(3). M is therefore a spherical triangle tiling the sphere, or, in the dihedral case, 
it is a lune (the situation discussed by Gromes). 

Making the special choice of conformal coupling in 1 + 2 dimensions means 
that the eigenvalues of the corresponding 'improved' Laplacian are perfect squares 
of integers, or half-odd integers, (for the unit sphere). By an appropriate selection 
of the action of the tiling group, T, one can arrange Dirichlet or Neumann conditions 
on the boundary of the fundamental domain, M.. A calculation, by separation of 
variables if you wish, e.g. [4] on the lune, gives the eigenvalues, A m , determined by 

u m = V\ m = a + d 1 m 1 + d 2 m 2 , m l7 m 2 = 0, 1, 2, . . . (4) 

where d 1 and d 2 are the integer degrees (not necessarily coprime) associated with 
the action of T, i.e. with the particular regular solid (3-polytope) . For example, 
for the lune of apex angle ir/q (q G Z), d 1 = q, d 2 = 1. The hemisphere corresponds 
to q = 1. 

The constant, a, takes the values 

a N = - and a D = d 1 + d 2 — a N 

for Neumann and Dirichlet conditions, [7]. 

One now sees the relation with Ehrhart polynomials. Trivially, counting eigen- 
values amounts to counting integer lattice points. It is no surprise that similar 
quantities appear. 
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Explicitly, consider the (special) rational triangle (2-polytope) 



V d = | (x 1 , x 2 ) € JR 2 : x i > 0, G^a^ + d 2 x 2 < 1 

with vertices (0,0), (l/d-^,0) and (0, l/d 2 ) where d 1 and d 2 are relatively prime, 
which is not an essential restriction. Then construct the integer dilation, lV d , and 
count the number, L(V d , I) of integer lattice points, (m 1 , m 2 ), inside (including the 
boundary as described previously). The result is given by (2) where now in fact c 1 
is a constant (because d 1 and d 2 are coprime). Hence, so far, 

- I 2 

L{VdJ) = 2d~d~ 2 +Cll + c ° {l) - 

The complete expression has been obtained by Beck and Robins, [2], and I 
reproduce the essential parts here, 

L(V d , I) = (l 2 + l(d 1 + d 2 + 1) + 1 (d 2 + d 2 + 3d,d 2 + 1) + <r(0) • (5) 

where a (I) is a periodic function involving Dedekind sums. The origin, Ehrhart [8], 
should also be consulted. 

To make a connection with eigenvalues, the essential counting restriction is 

d 1 m 1 + d 2 m 2 < I 

which, in terms of eigenvalues, (4), reads 

w m < I + a . 

so that the identification 

N(u) = L(V d ,u-a) (6) 

can be made but only when u — a is an integer. As a function of a real /, L(V d , I) 
will provide an interpolation and then the corresponding interpolated N(u) from 
(6) would give a smoothing of the exact N(u). To see how this compares with the 
smoothed expression (1), I put in the geometric values appropriate for a fundamental 
domain on the 2-sphere to get for (1), 

N(u) ~ ^ ± (di + d 2 - 1) w + J + . . . ) • (7) 
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The factor 2d 1 d 2 is the order of the reflective tiling group, V and the d± and d 2 need 
not be coprime. 

Choosing the upper (Neumann) sign for ease (a = 1/2), substitution of I = 
to — 1/2 in (5) does not yield (7). However, all is not lost. Plotting the exact lattice 
counting expression (5) at integer values of / produces a familiar staircase function. 
The interpolation provided by (5) for real / joins the upper part of the steps. The 
curve that joins the lower part is given by L(V d , I — I) and averaging these produces 
an interpolation (i.e. smoothing) that passes half way up the vertical rises. It is 
easily confirmed that this symmetrical combination 

±{L(V d ,u -1/2) + L(V d ,uj -3/2)) (8) 

reproduces the first two terms of (7), with the upper sign. This is the conclusion so 
far. 

Also one obtains the equality 

\M\ = 4n\V d \. 



3. Higher dimensions 

To show that this result possibly goes beyond mere numerical coincidence 
(which is unlikely anyway in view of the dependence on d 1 and d 2 ), I look at three 
dimensions more particularly. 

Generally, in d dimensions, the measures of the fundamental domain and its 
boundary are given by 

\M\ = ^\S d \ and \dM\ = b -^\S d ~ 1 \ 

in terms of the order, 2g, of the tiling group, T, and the number, 6 l5 of reflect- 
ing great hyperspheres (or reflecting <i-flats in the embedding (d + l)-dimensional 
Euclidean manifold). In terms of the degrees 

d d 
g = l[d l and b 1 = ^(d, - 1) + 1 

i=i i=i 

so that, in particular for d = 3, g = d x d 2 d 3 and b x = d 1 + d 2 + d 3 — 2. 
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Writing out the conventional asymptotic behaviour of the smoothed counting 
function in this particular situation gives 



where I have included the next term which involves b 2 , the number of elements of 
the group V that fix a (d — l)-flat but not a <i-flat. In three dimensions, 

b 2 = d 2 d 3 + d 3 d± + d±d 2 — d 1 — d 2 — d 3 . 

In three dimensions the Neumann constant, a N , equals 1 and the suggested 
symmetrical (midway) combination is, instead of (8), 

±(L(V d ,u-l) + L(V d ,u-2)). (10) 

The relevant polynomial terms have been determined by Beck et al, [9], and 
again I write them out 



l - (3(d 1 + d 2 + d 3 + d 2 d 3 + d^ + d x d 2 ) + d\ + d\ + df) + . . . ^ 



(11) 



where the omitted terms are periodic. 
Evaluating (8) yields the expression 



u 3 j u 2 (b l (b 1 -l)+b 2 + l)u 



2d x d 2 d 3 V 3 1 2 6 

the first two terms of which agree with the 'standard' expansion (9). The third term 
differs, but only very slightly, which proximity I cannot explain. One might expect 
that one should reproduce those polynomial coefficients which are constant, i.e. 
all those except the final one, proportional to w , although the higher heat-kernel 
coefficients depend on the propagation equation. 
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4. Degeneracies as denumerants 



I now wish to proceed to some exact relations, with the aim of making the 
previous conjecture more sensible. They are also more interesting. 

Averaging the counting function involves loss of information. I return to the 
exact function which determines, and is determined by, the (square root) eigenvalue 
spectrum, {u k }, k = 1, 2, . . ., coincidences implying degeneracy. 

The spectrum can also be specified by the distinct eigenlevels, say, / = 
0, 1,2, . . ., and the corresponding degeneracies, g(l). On the fundamental domain, 
leaving the eigenvalues in the form (4) expressed in terms of two integers (or more 
for higher spheres) one can often avoid explicitly constructing the degeneracies. 
However degeneracies have a useful role and appear in other contexts. 

For a given Neumann eigenlevel, u(l), the number of integer lattice points 
(m 1 ,m 2 ) satisfying 

d 1 m 1 + d 2 m 2 = I = — a , with a = - , (12) 
gives the degeneracy, i.e. 

g N (l; d) = |t ^(m l7 m 2 ) G Z+ : d x m x + d 2 m 2 = l^j . 

g^(l', d) counts the number of lattice points on the hypotenuse of the dilated ra- 
tional triangle lV d and is often referred to as a restricted partition or a denumerant, 
in an older terminology. Bell, [10], denotes the obvious extension to higher dimen- 
sions by D{1 | d) and Sylvester employs a variety of notations and terminology. His 
final one seems to be 

D(l | d) 



di, d 2 , . . . , d d , 
which I shorten to l/d ,. 4 

The study of denumerants is effectively the same as that of Ehrhart (quasi) 
polynomials as discussed later. Ehrhart himself, [8], actually spends most analytical 
time on them. 

The associated generating function is well known, going back at least to Euler, 

[11], 

CXD ^ 

h N (a; d ) = J2 9n(1; d) ° l = {1 _ a d l){1 _ a d 2) • ( i3 ) 



4 I retain the geometric terminology of 'degrees' for the denominator of the denumerant. Sylvester 
refers to them as the 'components' of the denumerant, among other things. Another possibility 
would be 'parameters'. 



6 



Of course, as a lattice statement, this holds for any integers d x and d 2 . For 
the S 2 /r spectral problem, d 1 and d 2 take only certain values and the generating 
function (13) is just the Molien series for a finite reflection group action, e.g. Stanley, 
[12], Meyer [13] and was explored by Laporte, [14], from a physical, mode point of 
view. The extension to higher dimensions is immediate. See Briining and Heintze, 
[6], Berard and Besson, [5]. 

Dirichlet conditions correspond to a different value for a in (12) and amount 
to a shift in I by a constant d where 

d = d 1 + d 2 -l, 

and is the number of reflecting hyperplanes in the embedding space. 
Then the Dirichlet generating function is 

h D (a; d) = g 9o(l; d) = ^—^—^ , (14) 

g D (l + d ;d) = g N (l;d), I = 0,1,2,.... (15) 
As an important example, for the lune of angle ir/q, q G Z, 

M(7;9 - 1)= (l-^Kl-a) 

M0r;g ' 1)= (l-^Kl-a) 
=h N (a;q, 1) - 



or 



1 -a 

A classic roots of unity, or trigonometric, calculation gives the explicit values 
for the degeneracies 

g N (l;q,l) = [l/q\ + l= \l/q\ , g D (l;q,l) = [l/q\, 

the solutions for the corresponding denumerants, (12), e.g. Sylvester [15]. 

Combining the Dirichlet and Neumann spectra (which amounts to adding spec- 
tral quantities such as the generating functions) gives the spectrum for a doubled 
fundamental domain on which the modes are periodic under the pure rotational 
part of T. Thus for the periodic lune (of angle 2n/q) we get the standard formula, 

1 + a q 

h(a; q, 1) =h N (a; q, 1) + h D (a; q, 1) = — - 

(1 - ai)(l - a) 

(16) 

= £(2Lf/gJ + iy, 

t=o 
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which exhibits the usual degeneracies sometimes obtained from group characters, 
e.g. [7], Harmer, [16]. 

For the finite number of periodic, uniform sphere tilings, the generating func- 
tions can be obtained by combining cyclic expressions corresponding to a geometric 
decomposition of the action of the subgroup of SO(3) in terms of axes and orders, 
[13,17,7]. Using the orbit-stabiliser theorem, this leads to the elegant result, [18], 



h(a;d) = ±(^2h(a;q,l)-h(<T-,l,l) y ) 



where the sum is over all axes of order q. As an example, the tetrahedral group has 
d x = 3, d 2 = 4 and q = 2, 3, 3. Simple algebra yields, using just knowledge of the 
qs, 

h(a;3,4) 



1 + 



(l-a 3 )(l-a 4 ) ' 

from which the values of d 1 and d 2 could be read off, if they weren't known. 

Using the cyclic degeneracies, (16), one derives the explicit expression for the 
rotational tetrahedral degeneracies, 



s(/;3,4) 



+ 2 



6 



+ 1- 



+ 1-/ 
- 2 



(17) 



For comparison, I reproduce the corresponding formulae for the octahedral and 
icosahedral tilings, [18], 



l 




/ 




/ 




+ 




+ 




2 


3 


4 



/ 

12 



+ 1 



+ 1 

l_ 

3 



(18) 



and 



</(Z;6,10) = 



/ 




/ 




/ 




+ 




+ 




2 


3 


5 



/ 



= 30 +1 



+ 1-/ 

/ 
3 



(19) 



Returning now to strictly lattice considerations, the solution of (13) for the 
lattice quantity, g N (t;d), with any coprime d x and d 2 , is given by Popoviciu's 
theorem, 



/ 

d x d^ 



+ 1 - 



d±H 



d<) 



d 2 H 
di 



(20) 
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where d 2 1 is the mod(<i 1 ) inverse of d 2 and d^ 1 that of d 1: mod^)- 

I wish to apply this to the tetrahedral case to check (17). I note that 3 _1 
3 and 4 _1 | = 1 so that (20) reads 



mod4 



SK(i;M ) = l + 1 -{!}-{! 



+ 1-/ 



3/ 




/ 




+ 






3_ 



(21) 



This is shown to be consistent with the rotational tetrahedral degeneracies, 
(17), by constructing the combination 

^(/;3,4) + ^(/-6;3,4) = i + ^-{| + i}-{|}-2{l 

according to the relation (15). Equivalence results in view of the equality, 

3/ li r 3/i r n 1 

essentially just a Hermite identity. 

This result means that one has explicit, separate formulae for the Neumann 
and Dirichlet mode numbers, discussed case by case by Laporte, [14]. I now check 
the other tilings, 

For this I require an extension of Popoviciu's theorem which is, (see Beck and 
Robins, [2] ex.1. 28) that, if gcd(d 1 ,d 2 ) = e, then 

Sjv ( ild ) = Jf r _(M\_/M\ + 1 , if e]l 



d x d 2 { d 2 J {d 1 ) (22) 
= , otherwise . 

where S 1 d 1 /e = 1 mod (d 2 /e) and S 2 d 2 /e — 1 mod (d 1 /e). 

In the octahedral and icosahedral cases, the greatest common divisor is 2 and 
I must be even for non-zero g N in (22). In the octahedral case, the reduced mod 
inverses are 8 1 = 2 and 5 2 = 1 so that 

<7iv(Z;4,6)=-^ + l-(H-(U, / even (23) 



12 [4 J L3^ 

which agrees with (18) because [l/'i] is zero for even / 
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The Dirichlet modes are obtained by applying (15) to give 



<? D (/;4,6) = 




again agreeing with (18), for I odd. 

A similar calculation holds for the icosahedral case where 5 1 = 2, 5 2 = 2. 
Substitution into (22) immediately yields equality with the icosohedral degeneracies, 
(19). 

In these two cases the rotational formula gives both the N and D values because 
these are associated with even and odd / respectively, whilst they are mixed up 
for the tetrahedron. A similar cross linking complication arises when computing 
spectral quantities for twisted vector bundles over homogeneous factors of the three- 
sphere, [19]. 

For the tetrahedron, because the degrees are coprime, there exists a Frobenius 
number which equals d l d 2 — d l —d 2 = 5, meaning that there is at least one A-mode 
for every integer / greater then 5. This is not true for the octahedron or icosahedron 
as the degrees are both even. 

5. Sylvester waves 

It is seen that all the relevant quantities take the standard form of a polyno- 
mial in / plus a periodic term, in agreement with an ancient partition theorem of 
Sylvester's, first stated in [15] to the effect that the denumerant takes the form, 

± = P(l) + U, (25) 

where P(l) is a polynomial in / of degree d — 1 (I work in d dimensions) and U, the 
'undulant' part, is periodic and contains roots of unity. 

Because of this periodicity, Sylvester refers to the parts of Z/d, in (25) as 'waves', 
depending on all roots of unity (most waves vanish by natural selection). 

For example, in the rotational tetrahedral degeneracy, (17), the periodic part 
has period 6 (the individual N and D parts having period 12). The leading term of 
the polynomial part is classic, [20]. 
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Sylvester, [15,21], gives the rule for calculating the polynomial part, P, as the 
coefficient of 1/t in the expansion of the quantity 



e it 



which is clearly a polynomial in /. 

By definition, this is directly expressed as a generalised Bernoulli function, 



(27) 



B M (-l\d) 



which is a polynomial in / whose coefficients are generalised Bernoulli numbers (up 
to a factor) . 

This expression for the polynomial part of the denumerant has been noted by 
Rubinstein and Fel, [22] who also point out the equivalence with a direct evaluation 
by Beck, Gessel and Komatsu, [23] which seems more or less the same as Sylvester's, 
[21]. 

The representation (27) is, in some ways, only cosmetically compact. It has 
value when general properties are being investigated and recursion relations etc. can 
be brought in. These can also be used to find explicit forms for particular cases, 
as discussed by Norlund, [24] and the Bernoulli polynomials can be computed in 
various ways. However, I prefer here to use the relation, employed in [25], with 
Todd polynomials, T n , [26], 



BW(x I d) = (-1Y r! $>l) s ^T;.^, . . . , a r _ s ) , r < d , (28) 



where the a s are symmetric functions of the degrees d i (i = 1, . . . , d). This form 
has the advantage of being valid for all d. That is, the functional form in the a i of 
the right-hand side does not depend on d. Accordingly, the expansion reads, 
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(d-iy-nAtS k r d - k ~ l[ l 
i d ~ 1 i k 

x ±z % i — n 



(29) 



k=0 

2 / 7d-2 7d-3 



n^A(^-l)! T ° + (rf-2)! Tl+ (rf-3) 



;T + 



Known forms for the Todd polynomials make the coefficients explicit, which 
allows checking, but the real advantage of using this structure is that the Todd 
polynomials are independent of d and one sees that the expression for dimension d 
is obtained by integrating that for (d— 1), the only 'new' term being the constant. 
I now present another (equivalent) version of this expansion. 

From an algebraic viewpoint a more advantageous variable is 7 = I + ^ Y2i 
Sylvester, [27] , calls this the augmented argument and also extends / to a real number 
so allowing the use of analysis, as, in fact, I have just done. 

A basic property of the Bernoulli functions, [24,28], shows that, expressed as a 
function of I, the polynomial part has the parity property, 

P d {-l) = {-^) d - 1 P d {l)- (30) 

Furthermore, the generating function equation takes on a more symmetrical 
aspect in terms of I, 

°° I - 1 

S ° l = TL (c- d ^ - <*'*) ' 



1=0 



which shows that, as a function of I, the full denumerant also obeys the parity 
property (Ehrhart reciprocity). It follows that the undulant part satisfies it also. 
This was proved by Ehrhart geometrically who thence derived the parity behaviour 
of the polynomial part. It seems that Sylvester, [27], had derived full reciprocity 
somewhat earlier and would also have had the reciprocity of Dedekind sums, if he 
had defined them. 

It is instructive to make the expansion (29) reflect the parity more easily. It 
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can be rewritten (using [24] p. 162, (18)), 

p ^ ( d - 1 )W. Ba(7+ ^ /2|d) 

in terms of the computable constants, [d] which vanish for odd z/. (This is the 
statement of reciprocity for this term). 

Again, the constants D^j [d] are essentially the polynomials, [26], and Hirze- 
bruch helpfully gives this time, the relation, 

D^[d] = ^A v {p 11 ...,p v ), 2u<d, (32) 

and the I expansion is, 

which terminates. 

Hirzebruch lists a few of the A polynomials as functions of the indeterminates, 
p i5 which are, here, the elementary symmetric functions of the squares of the degrees, 
[26], §1.3 which, to remind, are d in number. Thus, 

M = -\pu A 2 = ^(7p 2 1 -4p 2 ), A 3 = -^^(l6p s -Up 2Pl + 31pl). (34) 
I also reproduce some values of the B and D constants as given in [24], p. 167. 

nW_i n(»)__I, n(«)_l, + -V tf 2 d 2 ^ 

i<j 

where the s p are the sums of the p-th. powers of the degrees, d^ These expressions 
are easily shown to agree with those obtained via the Todd or the A-polynomials. 
Furthermore, they hold for all n and the conclusion is that one can relax the con- 
ditions in (28) and (32). Thus the range of 2u in (32) can be extended beyond d 
and, for example, from (34), one gets = 7df/15 agreeing with (35). 5 



5 More generally, the coefficient of the term in A v equals (2 2 "/(2i/)!) D 2v where the D 2v are 
listed in Table 4 in [28] and are the coefficients in the expansion of the characteristic series 
xj sinha;. 
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It should be apparent that there is no need to introduce generalised Bernoulli 
polynomials. It is possible to move to the final formulae in (29) and (33) directly in 
terms of multiplicative sequences. In fact this is what Sylvester does in his explicit 
determinations. 

The polynomial part, P is that Sylvester wave corresponding to the root of 
unity, 1. He denotes it by W 1 , and in his series of papers on partitions explains, 
several times, a systematic computational scheme for its algebraic determination, 
the most detailed being [27]. There he obtains, by rapid direct calculation, an 
expansion identical to (33), except that the polynomial coefficients are essentially 
all the homogenous products (of a given order) expressed in terms of functions of 
the sums of even powers of the degrees. I review the calculational details. A good 
introduction to Sylvester's theory can be found in the textbook by Netto, [29]. 

Instead of (26), one now has the more symmetrical form, 

e lt e ~it 
Wi = co i — j = co i — j . (36) 

nil (e d ^ 2 - e-W) -'ni^sinhi^ 

Expansion of the numerator yields the coefficient of the power / n /n!, as 
-^^sinhl^ Udi 



1 -r-r X 



Udi ^-i-llsinh^ 

where I have set x i = td i j c l and have encountered a multiplicative sequence, K, 
with characteristic function Q(x) = xj sinhx so that 

Q(x 1 )Q(x 2 ) ... = 1 + K 2 (s 1 ,s 2 )+ K A {s^ s 2 , s 3 , s 4 ) + . . . (38) 

where s k is the sum of the /cth powers of the x^ It is more conventional to use the 
symmetrical products, but this choice is more natural, as will be seen. The K m 
are homogeneous of degree m in the x i and so, in view of the definition of the x i , 
the right-hand side of (38) is a power series in t 2 . According to (37) the required 
coefficient is simply K d _ x _ n . To actually compute this, it is necessary to expand the 
particular characteristic function as a power series. Although everything is classical, 
I will do this ab initio from the old product 

x °° 1 

Q( x ) = = TT t, — 9/ ? o\ ■ ( 39 ) 

smhx Ai - (l + £ 2 /n 2 7r 2 ) 

n=l V ' / 



-n 
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Therefore 



\ogQ(x) = - ^ log (1 + x 2 /n 2 n 2 ) 

n=l (40) 



TV 2 " ' 2tt 4 "' v ; kn 2k 



The expansion of log Q rather than Q is the tactical point here and the con- 
struction of the multiplicative sequence, (38), follows on exponentiation, 



O O -expf - C(2) s t 2 + C(4) s t 4 - + tlnm s t ™ ± 

i .„ , 



™1>( v,/ 2 • -T 2 t 4 ... • ! i) k -r k t 2k : ... 



(41) 



I have introduced Sylvester's notation, 



m s (-!) fc+i ( 2k)c( (2k -i))s - ( ~ 1)fe+1 



on the modern definition of the Bernoulli numbers. In some ways it is better to 
leave these, and similar, constants in terms of (^-function values. 

The final step in determining the multiplicative sequence is to expand the right- 
hand side of (41) as a power series in t 2 . This can be done elegantly by relating 
this to the generating function of homogeneous products which is, (e.g. Littlewood, 
[30]), i 

i-h lX +h 2 x 2 -... = Y[— — 

i 1 

where h r is the sum of all homogeneous products of degree r of the indeterminates, 
x { . Taking logs, expanding that on the right-hand side, performing the sum on % 
and then exponentiating, to get back to the left-hand side, gives 

1 — h x x + h 2 x 2 — . . . = exp ( — s x x -\ — s 2 x 2 — . . . ) (43) 

and the h r are to be considered as functions of the sums of powers, s q , of the x^ 
Comparison with (41) shows that 

Q 1 Q 2 . . . = 1 - H x t 2 + H 2 t 4 - H 3 t 6 + . . . + (-l) r H r t 2r + . . . (44) 

where the H r are the same functions of the r i as the h r are of the s i . (A separate 
symbol is therefore not really needed.) 

15 



Solving the recursions that follow from the logarithmic derivative of (43) gives 
Brioschis's determinantal expression for the H r , (Faa de Bruno, [31] §89, [30]), 



H r = — 
r! 



Sylvester's final answer is, 

T d-i 



T l 


-1 











T 2 


T l 


-2 








r 3 


^2 


T i 


-3 





T r 


T r-1 









w 1 = 



1 



r 



T d-3 



T d-5 



UidAid-iy" (d-3)l Hl+ (d-5)! 



:tfo + . . . 



(45) 



(46) 



which is the same as (33) expressed slightly differently (H = 1). I observe that 
Sylvester was led naturally to multiplicative sequences. 

Sylvester, [27], lists six explicit (numerical) polynomials, which are also given 
by Ehrhart, [8]. 



6. The second wave 

There is a wave for every root of unity, I discuss in this paper only that one 
corresponding to the root —1, denoted by W 2 , which, like W 1 , can be algebraically 
expressed. Sylvester, [21], provides the following algebraic statement: W 2 equals 
the coefficient of 1/t in the generating function, 



;-i)' c 



IJt 



(47) 



where the a i are the even degrees, and the • are the odd ones. 

In order to express this in a form similar to (27), I recall the expansion formula 
for the generalised Eulerian functions e.g. [24,28], 



n„xt 



to which I add, for handiness, the one for Bernoulli functions, already used in (27), 



lGU(e< 



V = 
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Constructing the product in (47) the condition for the removal of the exponen- 
tial is 

1 a 

I = x + y--^a 4 --^/? 3 , (a + /3 = d), (48) 

i=l j=l 

which allows two simplest (in terms of I) solutions, 
and | 

the first of which yields the relevant coefficient required by Sylvester's rule as, 

w > = mJ-t ma , I (" * 1 ) Bi ° >(f+ Eai/2 1 a) E ' !l ~ M ' (49) 



i/=0 

while the second form leads to, 

a-l 



W > = 2^(t% ma , g ( V) TD -- M ^ (? + E/V2 ' ' (60) 

The constants ^ n) [/3] = £^ n) (£/V 2 | /3) a re zero for odd i/, [24]. 

Either form shows that W 2 is a 2-per iodic function of /, through (— 1) , mul- 
tiplied by a polynomial in I. For an odd (even) number, a, of even degrees, the 
polynomial has even (odd) order. If there is only one even degree, this polynomial 
is a constant, while W 2 is zero if all degrees are odd. (No degree has a factor of 
two.) 

It is possible to check reciprocity for W 2 directly. Changing the sign of I in the 
summand gives a factor of (— l) u which, because of the restriction on the index of 
the E, equals (— l) a_1 which can be taken outside the sum. Also the sign (— 1)' 
changes to (-1)"'"E Q -E ?i = (-1) Z "E^ = (-l) /_/3 so that the total change is 
to give a factor of (— l)'~^~ a + 1 = (—i)'-^ 1 i n place of (— l) z , which is the required 
reciprocity. 

Just for comparison, I give the formula that follows from an '/-solution', 

x = i + Yl a i y = ^2Pj 

of (48), viz, 

W 2 = 2P{ J:% Uai E ^ ~ ^ ( l + So* I «) ^-.(D*; I (3) • (51) 
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This form is given in the interesting paper by Rubinstein and Fel, [22]. It is, 
perhaps, less natural, reciprocity being harder to spot. Rubinstein and Fel also 
write the other waves in a similar way more complicatedly using the generalised 
Eulerian polynomials defined by Carlitz, [32]. 

I now derive a much neater expression for W 2 - As might be expected, Sylvester 
provides a route to an explicit solution, which I describe more compactly. I therefore 
begin again and am required to compute the coefficient, 



W 2 = (-l)'co_i 



e it 



n^sinhic^ H^coak^Pjt 



Similar to the W 1 computation, expansion of the numerator gives the coefficient 
of the power (— l) l J n /nl as, 



CO 



-l 



n,-2siiihito i n,-2coshit J 9, 
1 L P. 1 1 (52) 



CO 



2/3 Ho, -i-llsinh^ycoshr?,' 

where ^ = toj/2 and rjj = tfij/2. The product of two multiplicative series now 
arises. One has been encountered previously, see (37), (38), and I refer to it as the 
untwisted series. I now cover the other, defined by the characteristic function, 

^^^( 1+ «^7^) ' (53) 

in like manner 6 , and call this the twisted series. Thus, 

log!%) = _ y > + k|V!) yl _„. + ^<j^m^ ± , 

in terms of the Riemann-Hurwitz ^-function, C fl (s,w). 

The similarity of this to (40) shows that the formal algebra is just the same, 
only the numerical coefficients differ, with B k replaced by (2 2k — l)B k . So I define, 

_ C fl (2M/2) _ (-l) fc+1 (2 2fc - 1) r54 , 
^ " ( 27r )2fc s 2fc - 2(2l)! 2k 2k ' 1 j 



6 The products (53) and (39) can be considered as examples of the Mittag-LefHer theorem which 
notion allows one to generalise the approach. 
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and the multiplicative sequence this time is (cf (44)), 

fi x Q 2 . . . = 1 - H^) t 2 + H 2 {q) t A - H z {q) t 6 + . . . + (-l) r H r ^) t 2r + . . . . (55) 

According to (52) a product of multiplicative sequences is required. This can 
be done from the series (44) and (55) but is best performed first at the exponential 
level using (41) and, 

n,n 2 ... = exp^-qt 2 + ^ 2 t 4 -... + (-l) fc ^t 2fc ±...^) . (56) 

Doing this, and then expanding I find, as before, 

Q 1 Q 2 ...n i 2 ... = l-H 1 (^ + r)t 2 + H 2 ^ + T)t /i -... + (-iyH r ^ + T)t 2r + ... . 

(57) 

The final answer for the second wave is then formally similar to that for the 
first as a terminating series, 

W * = ^ l W^—{r-^\ H ^ - T~^ H ^ + T—*\ H * + •••)' (58) 
2P U.j<Xi V(«- !) ! (a -3)! (a -5)! J 

where H i = ii/^T + <;). T i defined by (42) with s 2q the sums of powers of the even 
degrees and q i is defined by (54) with, this time, s 2q being the sums of powers of 
the odd degrees. I could not find this exact form in Sylvester's writings. 

From a formal point of view, I note that the functions associated with the 
twisted multiple sequence are the generalised Euler polynomials. 

From a spectral perspective, the expansion variable is the eigenvalue, u, rather 
than /, hence I require 

(-r 1 * (*-<* I d) = (-I)'" 1 £ ( d " 1 ) (-a,)'5&_.(a I d) . (59) 



So that the polynomial, P, is 

d-i 

fc=0 



with 



Pfc = (rf _ 1 _ A:)! Jfe! Ili rfi Sd - 1 - fc(a 1 d) 
2 c 

~ r((fc + i)/2) °(d-fc-i)/2 

in terms of the heat-kernel coefficients, C*, on S d /r computed in [7]. 
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7. Ehrhart quasipolynomials 



Regarding the Ehrhart polynomial, I note that this is what could be termed 
an 'accumulated degeneracy', defined by, 

G(l) = 5>(0, (60) 

l'=0 

generically expressed. It is commonplace in generating function circles to write the 
obvious recursion for G as 

H(a) = Y ^—h(a) (61) 

for H(a) = J2Zo G ( l ) ° l etc - 

Exactly the present situation of spectra on regular tesselations of the <i-sphere 
has already been considered in [33] (extended to p-forms) and some discussion of 
the asymptotics were there made, using generating functions. 

Looking back at say (13), the recursion division in (61) is tantamount to adding 
1 to list of degrees, d { , and increasing d by one. 7 The denumerant for this case is 
the Ehrhart quasipolynomial L(P d ,/), i.e. 

The strictly polynomial part, P, of the Ehrhart quasipolynomial is given by 
the coefficient of 1/t in, 

„it 



nil (l-e-^)(l-e-*) 
which has already been found in (27) as (with a = (d — l)/2), 



\d d 



(62) 



d! 

The relation with Todd polynomials gives 
B^\x | d, 1) = (-1)- r! ^(-l) s ^T r _ s (a l5 . . . , a r _ s ) , r < d, (63) 



si 

s=0 



^ This is seen by noting that adding lm 3 to the left-hand side of (12) is equivalent to adding the 
results for (12) with the right-hand side equalling I, I — 1, 1 — 2, . . . , in turn. See [27], [8], [2]. 



20 



where now the a s are symmetric functions of the degrees d i (i = 1, . . . , d) and 1. 

First, I wish to express the symmetric functions, a s , of i = 1, . . . , d and 1 
in terms those, cr s , of just the d i since these are the variables used by [2], although 
this is not essential. Trivially, or from the defining fundamental identities, 

{[(i + d t t) = J2a r r 

i=l o 
d d 



H(l + dMl+t)=J2 W rt r , 



1=1 



it follows that 



a s = a s + a s _ x . (64) 
For convenience I list a few resulting Todd polynomials obtained from the usual 



ones, 



r„ = l, T 1 = i(a 1 + 1), T 2 = ^(a 2 1+ a 2 +3a 1 + l) (65) 



and write down some Bernoulli functions from (28), 
B[ d+1 \(d-l)/2\d,l) = ^±T -T 1 



B { 2 d+1 \{d- l)/2 | d,l) = 2T 2 - (d- 1)T X + ^-J-^j T 

The first two members of the polynomial (62) now read 
(-1)" f ujd+ d(2 + a 1 -d) uid _ 1+ ^ 



d 



(-i) rf f ld | dK + i) ^ | 



(66) 



& el d, 

where, for the combinatorialists, and as a check, I have transformed the series to 
one in /. These forms should be compared with (5) and (11). 

Of course it is not necessary to go via u to /, one can write (no different to 
(62)), d 

P ^irp- Br,(i|d ' 1) 

1 d l k 



A— 
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which is really only a re-expression. Using the forms of the Todd polynomials, (65), 
published coefficients are readily confirmed in rapid fashion. 

However, rather than use these expressions, it is far easier to take advantage 
of Sylvester's explicit expansion, (46) and, making the adjustments to move to the 
Ehrhart polynomial, I get for its polynomial part, 

1 ft f~ 2 f~ 4 \ 



YlidAdl u (d-1)! 1 (d-2)! 
where the H i are the same functions as before, (45) , but of new rs defined by, 

obtained by adding 1 to the list of degrees. The s q are sums of powers of the d { . 

Particular cases are easily constructed and the general one machine coded. 

I can now return to the recipe given in section 2 and construct the combination 
of Ehrhart polynomials, 

1 r ., (-l) d ( H d(a-,-d+l) i , 
- [A M + A(u - 1)] = + + ■ ■ • 

which agrees with the first two terms of the asymptotic expansion, (9), of the 
smoothed counting function, since b l = a x — d + 1 and confirms the statement made 
in section 2, for the first two terms. 



8. Conclusion 



It has been shown that the first two terms in the orthodox asymptotic expansion 
of the smoothed counting function can be obtained from a symmetrical combination 
of two Ehrhart polynomials in any dimension. 

It is, perhaps, no surprise that the heat-kernel coefficients can be obtained 
from the Ehrhart polynomial as this encodes all the eigenvalue information. 

The technical evaluation of denumerants is greatly eased by the use of multi- 
plicative sequences. Indeed these arise naturally. A simple expression for the second 
wave involving the numbers of all homogeneous products. In a later paper I intend 
to look at waves beyond the second one. 
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